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ABSTRACT 


An  invariant  second-order  closure  program  for  compressible 
shear  flows  has  been  developed  at  A.R.A.P.  The  program  has 
been  applied  to  study  a  simulation  of  the  mixing  region  in  a 
HF  chemical  laser.  Considerable  insight  into  the  complex 
processes  occurring  in  laser  cavit.ies  has  been  achieved  by 
these  studies.  The  results  show  that  the  assumptions  used 
in  the  eddy  viscosity  and  the  turbulent  kinetic  energy  closure 
methods  are  not  completely  valid  for  nonequilibrium  flows  such 
as  in  the  chemical  laser  system. 

The  results  of  studies  of  other  shear  flows  with  the  same 
invariant  second-order  model  are  in  generally  good  agreement 
with  experimental  measurements,  but  suggest  the  need  for  im¬ 
provement  of  the  pressure  diffusion  model  for  high  speed  com¬ 
pressible  flows. 

A  "typical  eddy"  box  model  has  been  developed  for  compu¬ 
tation  of  higher-order  correlations  involving  species,  temper¬ 
atures  and  reaction  rate  fluctuations  for  closure  of  the  turbu¬ 
lent,  multi-species,  reacting  flow  equations. 
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INTRODUCTION 


: 

; 


i . 

A  second-order  closure  program  for  the  investigation  of 
compressible  chemically  reacting,  turbulent  flow  fields  is 
under  development  at  A.R.A.P.  For  the  past  several  years 
A.R.A.P.  has  teen  in  the  forefront  of  the  application  of 
second-order  modeling  techniques  to  the  calculation  of 
turbulent  shear  flows.  These  techniques  have  been  success¬ 
fully  applied  to  the  calculation  of  planetary  boundary  layers'*', 
free  shear  layers  ,  axisymmetric  wakes  ,  flat  plate  boundary 
layers  and  other  turbulent  flow  problems.  The  new  program 
will  extend  the  methods  developed  in  these  studies  to  com¬ 
pressible  flow  problems  involving  highly  exothermic  chemical 
reactions  (i.e.,  the  turbulent  mixing  and  reaction  in  a  HF 
chemical  laser).  The  mixing  in  chemical  lasers  takes  place 
under  conditions  of  large  local  heat  release.  The  resulting 
interaction  between  the  turbulence  and  the  heat  release  is 
an  important  feature  of  the  flow. 

During  the  past  year,  with  support  from  AFOSR  Contract 
No.  F44620-73-C-0027 ,  we  have  employed  the  tehcnique  of 
second-order  modeling  to  study  various  aspects  of  gas  flows 
in  chemical  lasers  and  the  effect  of  Reynolds  number  and  wall 
roughness  on  the  mixing  in  the  laser  cavity.  Other  problems 
of  particular  interest  to  the  Air  Force  Weapons  Laboratory 
have  also  been  investigated.  These  studies  were  summarized 
in  the  Interim  Report  of  February  1974-5.  In  the  past  four 
months  a  number  of  additional  studies  on  the  laser  flowfield 
have  been  completed  and  are  discussed  in  this  report.  The 
shear  flow  program  has  also  been  used  to  check  some  of  the 
third-order  models  by  comparison  of  program  predictions  to 
available  experimental  data  for  certain  extensively  studied 
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shear  flows  such  as  the  flat  plate  boundary  layer  and  the  free 
shear  layer.  These  studies  are  presented  in  this  progress 
report . 

The  development  of  the  reacting  flow  program  is  pro¬ 
ceeding  rapidly.  The  modeling  of  the  third-order  correlations 

such  as  k  c^  ,  k'c^  c^  which  appear  in  the  transport 

equations  for  the  second-order  correlations  employs  a  "typical 

eddy"  box  model.  Extensive  testing  of  the  model  has  been 

carried  out  with  very  good  results.  The  model  is  described 

in  Section  2.5  of  this  report.  The  TEN SR  program  takes  the 

basic  transport  equations  in  general  tenser  notation  and 

converts  and  expands  them  to  the  corresponding  equations  in 

the  desired  coordinate  system.  It  was  originally  written 

only  for  Cartesian  coordinates.  It  is  being  improved  and  extended 

to  generalized  coordinate  systems  Including  cylindrical  coordinates. 

A  brief  description  of  the  program  is  piesented  in  Section  z.4. 

2.  PROGRESS  REPORT 
2 . 1  Laser  Cavity  Flow  Studies 

The  simulation  of  the  mixing  processes  in  a  HF  chemical 
laser  cavity  using  cur  turbulent,  -ompressible ,  shear  flow 
program  showed  that  at  t.he  typical  operating  conditions 
( Re ~  2000  for  the  AFWL  HF  laser)  the  turbulence  decays  rapidly. 

Under  these  operating  conditions,  the  performance  of  the  laser 
would  be  strongly  affected  by  changes  in  flow  conditions.  The 
turbulence  in  the  region  of  mix  i  rig  and  chemical  reaction  is 
simply  the  decaying  turbulence  left  over  from  the  wall  boundary 
layers  and  the  recirculation  regions  at  the  nozzle  exit.  These 
studies  also  showed  that  Reynolds  numbers  of  the  order  of  lO'’ 
would  be  necessary  to  reach  significant  turbulence  levels  in 
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the  mixing  region  which  may  lead  to  an  improvement  in  the  laser 
performance , 


A  comparison  of  our  second-order  closure  model  with  the 


eddy  viscosity  and  the  kinetic  energy  approach  has  been  made 
for  the  laser  flow  problem.  Figure  1  shows  the  results  for 
the  computation  of  ar.  eddy  viscosity  "constant" 


\uv  /max 

_ _ aJ _ 

tu  \  (y  -  y^) 


using  the  results  for  the  Reynolds  stress  u  V  and  the  mean 
axial  velocity  u  obtained  from  the  second-order  closure 
program.  We  see  that  K  is  neither  a  constant  with  the 
downstream  distance  X  nor  with  the  Reynolds  number,  and 
the  use  of  a  simple  eddy  viscosity  in  the  mixing  region  of 
interest  for  a  chemical  laser  would  be  incorrect.  The  figure 
also  shows  the  values  for  K  oMainea  for  a  free  shear  layer 
for  which  we  see  that  after  a  short  initial  region  K  aces 
become  a  constant.  Therefore,  in  an  equi 1 i b r lum  flow  region, 
the  eddy  viscosity  concept  can  be  used  with  some  degree  of 
success,  but  for  nonequilibrium  flows  as  in  a  lase>.  cavity  a  more 
sophisticated  procedure  lik>-  sec<  nd-<  rder  insure  must  be  used. 

An  elementary  second-order  closure  procedure  is  the 
turbulent  kinetic  energy  closure  method  rs  developed,  for 
example,  by  Harsha  .  Figures  2a  and  2b  show  the  results  of 

our  computation  of  the  quantity  a,  -  for  tl  e 

1  q 

laser  flow  field  at  Re  =  2x10  '  and  ?xl(.  .  We  see  that 

can  be,  to  a  large  extent,  treated  as  u  constant  across 


KUHN 


Free  shear  layer 
value 


Fi  (ju  r'e  1 


Variation  of  Iho  ►•dfiy  viroor  1  t,,y  ' 
rnixitir,  rof.ion  ol'  a  rtl  laia.'j 


the  flow  field.  There  appears  to  be  some  consistent  variation 
in  the  edge  regions  of  the  flow,  but  there  the  reliability  of 
the  calculations  is  poor  as  both  u'v'  and  q  are  small  in 
these  regions.  The  results  show  some  effect  of  the  flow 
Reynolds  number  with  a^  -  0.20  a,t  Re  =  2x10“  and  a,  -  0.12- 
0.16  at  Re  =  ZxlO'’.  Harsha  recommends  a  value  of  0.12  from  com¬ 
parisons  with  different  equilibrium  shear  flows .  In  general, 
then,  a  kinetic  energy  closure  method  seems  to  be  considerably 
more  successful  than  an  eddy  viscosity  model  but  still  falls 
short  of  being  able  to  predict  the  results  obtained  by  a 
second-order  closure  procedure. 

Our  second-order  closure  program  provides  detailed 
information  on  turbulence  components  useful  for  studying  the 
behavior  of  the  cavity  flowfield,  for  example,  the  study  of 
laser  beam  degradation  due  to  T  7T  ’  fluctuations  presented  in 
an  earlier  report.  A  detariieu  3eoond-order  closure  program 
has  to  be  developed  in  order  to  •  rrectly  incorporate  the 
interaction  between  the  turbulence  and  the  chemistry  in  a  HF 
laser.  We  are  continuing  cur  development  of  such  a  program. 

2.2  Shear  Flow  Studies 

An  Invariant  second-order  closure  model  requires  the 
basic  assumption  that  third-O:  -ier  and  higher-order  correlations 
that  appear  in  the  transport  equations  for  the  turbulence 
stresses  can  be  modeled  in  terms  of  the  second-order  corre¬ 
lations  and  mean  properties  in  ..o  h  a  way  that  tire  model  remains 
invariant  with  respect  to  changes  in  the  flow  geometry.  Models 
have  been  developed  by  ns  for  various  third-order  correlations. 
These  models  are  being  tested  and  refined  by  comparison  of 
model  predictions  bo  reliable  experimental  measurements  in 
various  shear  flows.  A  compilation  of  the  models  being  used  in 
our  present  programs  is  presented  in  reference  2.  Some  of  the 


results  of  our  studies  on  the  free  shear  layer  and  the 
flat  plate  boundary  layer  are  briefly  discussed  below. 

Figure  3  compares  the  predictions  for  the  effect  of 
the  velocity  ratio  on  the  spreading  parameter  c  for  an 
incompressible  free  shear  layer  against  experimental  measure¬ 
ments  and  a  generally  accepted  empirical  curve.  The  results 
are  in  excellent  agreement.  Figures  4  and  5  presents  compar¬ 
isons  of  our  program  predictions  against  Coifs'  law  of  the 
wall  profile  and  correlation  of  skin  fricti'.  n  data  for  in¬ 
compressible  flat  plate  boundary  layers .  'he  results  again 
are  in  very  good  agreement,  fne  model  predictions  for  the 
detailed  turbulent  correlations  are  compared  to  Klebanoff’s 
measurements  in  Figure  6.  The  agreement  for  u  'v '  i.s  very 

good.  The  results  f or  y/ u  no  not  show  the  sharp  peak  in 
the  region  very  near  the  wall,  but  the  agreement  is  good  foi' 

v/u^*  and  77^  over  most  ci  tr.e  profile.  Our  ore  sent 


v  u  *  and  y/  V  over  most  cl  tr.e  profile.  Our  present 
model  predicts  virtually  the  came  values  fur  v'  ’  and 


while  the  measurements  snow  a  significant  difference. 

It  must  be  pointed  out  that  other  measurements  cf  y/ 
indicate  considerable  uncertainties  in  tee  \  dues  for  this 
term  but,  nevertheless,  our  in:  .r.p  ress.lt, le  rac.iel  needs  some 
modifications  for  this  aspect  1  the  experimental  data. 

The  present  compressible  meted  neglects  treasure  diffusion 
effects  and  a  recent  study'  seems  to  indicate  that  it  might 
become  important  even  at  moderate  Mach  number...  Some  of  the 
disagreement  between  the  model  predictions  and  the  experimental 
data  for  compressible  flows  presented  bed, w  could  be  due  to  the 
pressure  diffusion  term.  We  are  presently  modifying  the  progran 
to  incorporate  our  pressure  diffusion  model  and  to  determine  the 
effects  of  this  term. 


O  Present  study  (o"0  =  ll) 

Y/A  Experimental  data  compiled 
by  Birch  &  Eggers 


Figure  4.  Mean  velocity  profiles  in  the  incompressible  flat  plate  boundary  layer. 


Figure  7  shows  the  effect  of  Mach  number  on  the  spieading  rate 
for  a  free  shear  layer.  The  predicted  Mach  number  effect  is  smaller 
than  suggested  hv  the  experimental  measurements.  This  is  believed 
to  be  partly  due  to  the  difference  in  the  velocity  ratios  (numerical 
;  roblems  in  the  computer  program  do  not  allow  shear  layer  runs  for 
u^/u^  <  0.1)  and  partly  due  to  the  neglect  of  pressure  diffusion 

effects.  The  results  of  Oh'  support  this  conclusion. 

r—  " 

Figure  8  compares  the  prediction  for  a  c  for  a  compressible 
boundary  layer  against  very  recent  measurements  cy  Johnson  .  This 
result  also  suggests  that  further  model  improvements  are  necessary 
for  compressible  flows  to  predici  details  of  the  flow  The  present 
model  does  an  excellent  Job  of  predicting  tr.e  effect  of  Mach  number 
on  the  skin  friction  as  shown  in  Figure  9,  where  model  results  are 
compared  with  previous  experimental  and  analytical  studies. 

2 . 3  Reacting  Flow  Program 

"Typical  Eddy"  Box  Model 

The  set  of  second-order  eiosu”e  equation;'  governing  the  mixing 
and  reaction  of  a  compressible,  turbulent,  mulct  component  flow 
system  contain  a  number  of  new  and  diffi  it  third-order  correlations, 

such  as  k  ’c  Jca '  ,  k  ;c  T  ,  k  V.  c'  ,  et  .  These  terms  being 

dp  cx  u  ri 

correlations  of  scalar  quantities  are  somewhat,  i:  ore  difficult  to 
model  than  correlations  involving  e  or  terns.  During  the  past 
year,  with  additional  support  from  NASA,  we  have  developed  a  "typical 
eddy"  box  model  that  is  very  promising  for  the  modeling  of  these 
terms  and  the  closure  of  the  reacting  flow  equations. 


The  model  can  be  looked  at  in  two  ways.  First,  one  might 
consider  the  model  to  be  a  straightforward  way  of  expressing  the 
higher-order  correlations  of  scalar  fluctuations  that  occur  in  the 
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Figure  8.  Jistribution  of  the  normalized  velocity  fluctuations 
across  the  boundary  layer 
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equations  for  the  second-order  correlations  of  fluctuating  scalars 
in  terms  of  the  second-order  correlations  and  the  means  of  the 
scalars.  Second,  it  can  be  looked  at  as  a  way  of  generating  the 
most  important  features  of  any  distribution  function  for  the 
turbulence  quantities  that  affect  the  chemical  behavior  0'  he 
system. 

The  model  we  have  selected  is  the  following: 

We  define  the  most  typical  eddy  that  will  pass  a  point  in  the  re¬ 
acting  flow  at  any  time.  For  the  purpose  or  computing  scalar 
correlations,  the  extent  of  this  eddy  in  time  is  unimportant  and 
is  taken  equal  to  one.  Consider  a  chemical  reaction  in  which  two 
species,  say,  a  and  p,  react  to  form  a  third  species  7  . 

The  typical  eddy  of  extent  unity  in  time  has  the  following 
structure.  For  a  time  out  of  the  total  unit  time,  the  eddy 

contains  only  the  species  a  .  For  a  time  equal  to  ,  the 

eddy  contains  only  the  species  p  .  For  a  time  equal  to  e  , 
the  eddy  contains  only  the  species  7  .  For  times  equal  > 

,  and  ,  the  eddy  contains  |  a  and  j  p  ,  =>-  a  and 

?  7  ,  and  ^  P  and  ^7  by  mass,  respectively,  in  a  state  of 
molecular  mixedness.  In  addition,  for  a  time  equal  to 


1 


the  eddy  will  contain  a  ,  p  ,  and  7  uni formiy  mixed  on  the 
molecular  level  according  to  the  mean  mass  fractions  a  ,  B  , 

and  7  . 


A  typical  eddy  may  then  be  drawn  as  in  Figure  10a  insofar 
as  the  partitioning  of  time  in  the  eddy  into  various  species  is 
concerned.  Note,  since  the  total  time  of  the  eddy  is  unity,  that 
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■en** 


Thus,  there  are  only  five  unknowns  necessary  to  define  the  eddy. 

Five  equations  are  available  to  determine  these  unknowns . _ Inde- 

pendent  equations  are  available  for  a  ,  3  ,  a'0'  ,  a'y'  . 


and  3'y',  and,  therefore,  we  are  able  at  ail  times  to  define 
an  eddy  in  terms  of  only  the  mean  arid  second-order  correlations 
that  are  being  calculated.  This  model  allows  calculation  of 
any  higher-order  correlations  of  .ass  fractions  in  terms  of 
a  ,  3  ,  a' 3'  ,  a  y'  ,  and  3'y  . 

To  proceed  with  the  calculation  of  other  higher-order 
correlations  involving  the  s  alai  T  (temperature),  we  make 
the  following  additional  assumptions.  For  the  time  t.  when 
the  eddy  contains  only  a  tne  base  temperature  of  the  eddy 
is  T  ;  for  the  time  c  when  the  eddy  contains  all  0  , 
the  base  temperature  is  ;  for  the  time  e ...  when  the  eddy 

is  all  y  ,  the  base  temperature  is  T._  .  Fci  the  times 
when  the  eddy  consists  of  i  a  and  ‘  o  b.  mass,  the  base 
temperature  is  taken  to  be  the  appropriate Ly  weighted  mean  of 
Ta  and  T^  .  We  do  likewise  for  cells  containing  \  a  and 

i  y  and  for  those  containing  -■  [3  and  ,  h,v  mass. 

We  also  do  likewise  for  the  cell  containing  a  ,  3  ,  and  y  . 
Since  the  model  must  admit  temperature  flu  tuarions  whose  dis¬ 
tribution  functions  are  almost  Gaussian  when  only  one  species 
is  present,  we  will  assume  that  each  ce'l  ; pends  one-half  its 
time  at  a  temperature  AT  above  its  rase  temperature.  This 
temperature  model  may  be  drawn  as  shown  in  Figure  10b.  We 
see  that  there  are  four  unknowns  •  this  temperature  model: 

Ta  *  Tj3  ,  T  ,  and  AT.  There  are  four  equations  to  solve 
these  quantities,  for  the  eddy  mm.:,  h  .  •.  th  >  average 
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temperature  T  ,  and  the  correlations  T'  .  a'T  ,  and  &'t' 
that  are  available  from  the  equations  that  are  carried  along  to 
define  the  flow. 

The  concentration  mode]  and  the  temperature  model  just 
descrioed  allow  one  to  compute  any  higher-order  correlation 
of  scalar  fluctuations  involving  concentrations  or  temperatures. 

The  second- and  higher-order  correlations  involving  the  fluctuations 
of  reaction  rate  constant  k'  can  also  be  calculated  from  the 
model  just  given.  TK.Vs  is  so  because  for  eaen  half  cell  shown 
in  Figure  1C,  it  is  possible  to  calculate  the  instantaneous 
reaction  rate  k  so  that  the  reaction  rate  distribution  may 
be  plotted  as  in  Figure  10c .  Fror:  this  figure  it  is  possible 
to  calculate  k'a'  ,  k' }  k'y'  ,  or  any  other  higher-order 
correlation  of  k'  with  fluctuations  in  concentration  or 
temperature . 

The  above  simple  model  has  01 onst  u  ted  after  studying 
a  number  of  other  models  and  ■  ■  ]ct  i  the  results  of  the  model 
predictions  to  exact  solutions  of  the  reaction  equations.  Details 
of  these  studies  are  available  li  aefei  eij.'e  9*  The  final  model 
selected  is  believed  to  be  capable  of  providing  the  proper  closure 
of  the  equations  and  to  make  second  -  •,?  dei  closure  bomputation  of 
exothermic  reacting  flows  possible  for  the  : ;  si  time. 

2 .  U-  TENSR  Program 

In  1970  when  A.R.A.P.  firs  app  cached  the  task  of  numeri¬ 
cally  solving  the  modeled  equations  for  compressible  turbulent 
flow,  it  was  decided  that  portions  of  the  process  of  producing 
the  required  program  should  be  automated.  In  this  way,  much 
tedious  algebra  and  detailed  coding  could  be  handled  automat¬ 
ically  and  a  major  source  of  errors  eliminated. 


Two  program  systems  were  developed  to  accomplish  this  task. 

Tne  first  program,  TENSR,  accepted  differential  equations  in 

tensor  notation  as  input.  The  program  expanded  the  terms  accord- 

ing  to  the  summation  convention  (Am  Bm  -*■  A~  B-^  +  A  +  A^  B^)  , 

and  reproduced  separate  versions  of  the  equations  for  each 

specified  value  of  the  free  index  or  indices  (C  .  =  D.  -*  C  ,  = 

Df,  C  2  =  I^,  C  ^  =  D^)  .  In  addition,  TENSE  did  extensive 

checking  of  the  validity  of  the  inputs,  allowed  the  substitution 

of  symbols  (U1.  =  0  —  U1,  +  UL\  f  \j\  =  0-*U,+V;>  +  W,=  0), 

,i  ,u  ,u  ,u 

and  evaluated  individual  components  of  the  Kronecker  delta  and  the 
metric  tensor.  The  output  of  TENSR  was  edited  by  hand  to  remove 
terms  not  applicable  to  the  problem  at  hand,  i.e.,  terms  that  were 
negligible  under  the  boundary-layer  assumptions. 

The  resulting  equations  were  inputs  to  the  second  program 
system,  DIFFR,  which  constructed  FORTRAN  statements  for  com¬ 
puting  the  coefficients  of  the  finite-difference  equations 
corresponding  to  the  input  differential  equations. 

These  programs  were  very  successful.  However,  there  was 
a  major  restriction.  Covariant  differentiation  was  taken  to 
be  the  same  as  partial  differentiation.  In  other  words,  the 
Christoffel  symbols  were  taken  to  be  zero.  The  results  were 
valid  only  for  systems  for  which  all  the  elements  of  the  metric 
tensor  are  constants.  Thus,  the  only  applications  were  for 
Cartesian  coordinates. 

For  the  present  effort  it  became  necessary  to  solve  even 
more  complicated  equations  in  cylindrical  coordinates.  There¬ 
fore,  the  TENSR  system  has  been  revised  so  that  it  automatically 
produces  the  Christoffel  symbols  associated  with  covariant  dif¬ 
ferentiation.  Orthogonality  was  not  assumed,  so  that  it  is 
possible  for  example,  to  introduce  streamline  coordinates  by 
defining  the  metric  appropriately. 
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In  making  the  revisions,  other  changes  are  being  made  in 
TENSR  as  experience  with  it  has  shown  these  changes  to  be 
desirable.  These  include  changing  the  verification  section 
to  make  the  verification  routine  more  efficient.  Another 
change  will  make  it  possible  to  automatically  accomplish  at 
least  some  of  the  editing  to  remove  negligible  terms  (mentioned 
above),  if  TENSR  is  so  instructed. 

No  major  changes  have  been  made  in  DIFFR,  nor  are  any 
contemplated.  However,  a  successor  program  has  been  written 
which  produces  more  efficient  FORTRAN  coding  from  that  produced 
by  DIFFR. 

3 .  CONCLUSION0 

Considerable  insight  into  the  complex  processes  occurring 
in  lasei  cavities  has  been  achieved  through  the  use  of  an  in¬ 
variant  second-order  closure  model  of  the  turbulent  flow 
equations.  The  assumptions  used  in  the  eddy  viscosity  method 
and  the  kinetic  energy  method  are  not  completely  valid  for 
nonequilJ brlum  flows  such  as  in  the  chemical  laser  system. 

Studies  of  different  incompressible  shear  flows  with  the 
same  basic  secord-order  model  demonstrates  that  the  general 
features  of  such  flows  are  described  very  well  by  an  invariant 
model.  It  is  believed  that  with  a  little  more  modeling  effort 
for  compressible  flows  (especially  for  the  pressure  diffusion 
terms)  all  the  mpcrtant  details  of  compressible  flows  can  also 
be  adequately  calculated  by  our  second -order  closure  scheme 
without  abandoning  the  notion  of  invariance. 

Significant  progress  has  been  made  towards  a  reacting 
flow  program  with  tire  development  of  the  "typical  eddy"  box 
model  for  computation  of  correlations  involving  species, 
temperature  and  reaction  rate  fluctuations. 
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